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Abstract 

We define a non-iterative transformation method for Blasius equation 
with moving wall or surface gasification. The defined method allows us to 
deal with classes of problems in boundary layer theory that, depending on a 
parameter, admit multiple or no solutions. This approach is particularly con¬ 
venient when the main interest is on the behaviour of the considered models 
with respect to the involved parameter. The obtained numerical results are 
found to be in good agreement with those available in literature. 
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1 Introduction. 


The problem of determining the steady two-dimensional motion of a fluid past a 
flat plate placed edge-ways to the stream was formulated in general terms, accord¬ 
ing to the boundary layer theory, by Prandtl 0]|, and was investigated in detail by 
Blasius El. The engineering interest was to calculate the shear at the plate (skin 
friction), which leads to the determination of the viscous drag on the plate, see for 
instance Schlichting If3l). 

The celebrated Blasius problem is given by 
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where / and 77 are suitable similarity variables and in literature we can find either 
P — 1/2 or P — 1. This is a boundary value problem (BVP) defined on the semi¬ 
infinite interval [0,°°). According to Weyl 0, the unique solution of (Q]) has 
a positive second order derivative, which is monotone decreasing on [ 0 ,°°) and 
approaches to zero as 77 goes to infinity. The governing differential equation and 
the two boundary conditions at 77 = 0 in (|T|) are invariant with respect to the scaling 
group of transformations 


fl* = A - “, f = l a f (2) 

where a is a nonzero constant: Topfer used a — 1/3, see 01, but we have always 
put a — 1 in order to simplify the analysis. The mentioned invariance property 
has both analytical and numerical interest. From a numerical viewpoint a non¬ 
iterative transformation method (ITM) reducing the solution of (Q]) to the solution 
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of a related initial value problem (IVP) was defined by Topfer 0. Owing to that 
transformation, a simple existence and uniqueness Theorem was given by J. Serrin 
[f6ll as reported by Meyer |71 pp. 104-105] or Hastings and McLeod | 8 ] pp. 151- 
153]. Let us note here that the mentioned invariance property is essential to the 
error analysis of the truncated boundary solution due to Rubel 0, see Fazio IflOil . 

Our main interest here is to extend Topfer’s method to classes of problems in 
boundary layer theory involving a physical parameter. This kind of extension was 
considered first by Na [Hill , see also Na Ifl2l Chapters 8-9]. The application of a 
non-ITM to the Blasius equation with slip boundary condition, arising within the 
study of gas and liquid flows at the micro-scale regime [13] fl4l . was considered 
already in [|T3Tl . Here we define a non-ITM for Blasius equation with moving 
wall considered by Ishak et al. llT6il or surface gasification studied by Emmons 
liTTI and recently by Lu and Law I1T81 . In particular, we find a way to solve non- 
iteratively the Sakiadis problem [ 19 J( 20]. For the solution of the Sakiadis problem 
by an ITM see Fazio ll2D . The defined method allows us to deal with classes of 
problems in boundary layer theory that, depending on a parameter, admit multiple 
or no solutions. This approach is particularly convenient when the main interest is 
on the behaviour of the considered models with respect to the involved parameter. 


2 Moving wall 


According to Ishak et al. |[T6l the differential problem governing a moving wall, 
with suitable boundary conditions, is given by 


d*£ , I 

dr \ 3 2^ dr] 2 
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where P is a non-dimensional parameter given by the ratio of the wall to the flow 
velocities. Blasius problem ([[]) is recovered from ([3]) by setting P — 0. 
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2.1 The non-ITM 


The applicability of a non-ITM to the Blasius problem © is a consequence of 
both: the invariance of the governing differential equation and the two boundary 
conditions at 77 = 0 , and the non invariance of the asymptotic boundary condition 
under the scaling transformation ([2]). In order to apply a non-ITM to the BVP © 
we consider P as a parameter involved in the scaling invariance, i.e., we define the 
extended scaling group 

r = A/, T]* = X~ x r\, P* = A 2 P. (4) 


Let us notice that, due to the given second boundary condition at 77 = 0 and the 
asymptotic boundary condition in ©, P has to be transformed under the scaling 
group © with the same law of ^(rj). By setting a value of P*, we can integrate 
the Blasius equation in © written in the star variables on [0,T]^j, where 77 ^ is a 
suitable truncated boundary, with initial conditions 

/*(0) = 0, $£(0)=r, g(0) = ± l, (5) 

in order to compute an approximation ^*( 77 ^) for ^*(°°) and the corresponding 
value of A according to the equation 


A = 


df* 

dr]* 


«)+P* 


1/2 


( 6 ) 


Once the value of A has been computed, by equation ©, we can find the missed 
initial conditions 


d£ 

dr\ 


( 0 ) 


A 2 P* 


dr] 2 


(0) = A 


-3 


d 2 P 

dr]* 2 


( 0 ). 


(7) 


Moreover, the numerical solution of the original BVP © can be computed by 
rescaling the solution of the IVP In this way we get the solution of a given BVP 
by solving a related IVP 

We remark here that the plus (for P < 0.5) or minus (when P > 0.5) sign 
must be used for the second derivative in ©. Moreover, the computation of a 
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value at infinity is unsuitable from a numerical viewpoint and therefore we use 
a truncated boundary 77 ^ instead of infinity. For the application of the method 
defined above, depending on the behaviour of the numerical solution, we used 

= 10 or 77 * = 15. 

In tableQ]we list sample numerical results obtained by the non-ITM for several 
values of P*. Here the D notation indicates that these results were computed in 
double precision. As mentioned before, the case P* — P = 0 is the Blasius problem 
©. In this case our non-ITM becomes the original method defined by Topfer |5J. 
For the Blasius problem, the obtained skin friction coefficient is in good agreement 
with the values available in literature, see for instance the value 0.332057336215 
computed by Fazio Il22l or the value 0.33205733621519630, believed to be correct 
to all the sixteen decimal places, reported by Boyd [|23l . The values shown in the 
last line of table |Tj are related to the Sakiadis problem [|T9l 1201 and were found 
by a few trial and miss attempts. For this problem, the obtained skin friction 
coefficient is in good agreement with the values reported by other authors, e.g. 
-0.44375 Sakiadis JED, -0.4438 by Ishak et al. [[13, -0.44374733 by Cortell 
Il24l or —0.443806 by Fazio lf2T1 . 

Figure Q] shows the solution of the Sakiadis problem, describing the behaviour 
of a boundary layer flow due to a moving flat surface immersed in an otherwise 
quiescent fluid, corresponding to P — 1. Actually, this is a case of practical interest 
if we are considering the plate as an idealization of an airplane wing. Let us notice 
here that by rescaling we get 77 /j, < 77 ^. 

,2 f 

In figure [2] we plot 4^4(0) versus P. From this figure we realize that our 
problem has an unique solution for P > 0, while dual solutions exists for P < 0. 
The values of the skin friction coefficient are positive for P <1/2 and they become 
negative when P > 1/2. From a physical viewpoint, this means that a drag force 
is exerted by the flow on the plate in the first case, while in the second the force 
is of opposite type. Figure [2] also shows that the solutions exist until a critical, 
negative, value of the parameter, say P c , is reached. From the data in table Q] we 
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d 2 f* 
%( 0 ) 
dr[*~ 

p* 

dr H 

dri * v ’ 

> 

P 

1 

-500. 

1.55D04 

5.46D-07 

-0.033393 


- 100 . 

2.34D03 

9.42D-06 

-0.044591 


-5. 

36.325698 

0.005704 

-0.159613 


-1.5 

4.368544 

0.205830 

-0.522913 


-1.25 

3.529165 

0.290627 

-0.548447 


- 1 . 

2.917762 

0.376537 

-0.521441 


-0.75 

2.503099 

0.430814 

-0.427814 


-0.5 

2.250439 

0.431797 

-0.285643 


0 . 

2.085393 

0.332061 

0 . 


1 . 

2.440648 

0.156689 

0.290643 


5. 

5.771518 

0.028287 

0.464187 


100 . 

1 . 00 D 02 

3.53D-04 

0.499557 


500. 

5.00Z)02 

3.16D-05 

0.499960 

-1 

100 . 

99.822681 

-3.54D-04 

0.500444 


10 . 

9.433763 

-0.011673 

0.514568 


5. 

4.182424 

-0.035939 

0.544519 


2 . 

0.528464 

-0.248722 

0.790994 


1.719 

-4.73D-05 

-0.443715 

1.000027 


Table 1: Moving wall boundary condition: non-ITM numerical results. 


get P c —0.548447. This value is in good agreement with the value P c = —0.5483 
computed by Ishak et al. [H6ll using an iterative method: the second order Keller’s 
Box finite difference method, see Keller 11251 l26ll . The boundary layer separate 
from the surface at P — P c , and, therefore, for P < P, the Navier-Stoker equations 
has to be solved because the hypotheses of boundary layer theory felt down. We 
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Figure 1: Numerical results of the non-ITM. Top frame: solution of the IVP; 
bottom frame: solution of the Sakiadis problem found after rescaling. 

have separation for a positive value of the skin friction coefficient and not at the 
point where this coefficient vanish as in the classical boundary layer theory. The 
zero value of the skin friction coefficient when P — 1/2 corresponds to equal 
velocity of the plate and the free stream and does not mark separation. 

3 Surface gasification boundary condition 

In this section we consider a surface gasification flow that, with and without sub¬ 
sequent gas-phase flame-sheet burning, was first formulated and analysed by Em- 
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Figure 2: Blasius problem with moving wall boundary conditions. Non-ITM: 
0 ( 0 ) versus P. 


mons da. For such a flow we have to consider the variant of the celebrated 
Blasius problem 


d*f_ 

dr \ 3 
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/(°) = 0 ) , 


df_ 

drj 


( 0 )= 0 , 


df_ 

dr\ 


(v) 1 


( 8 ) 


as 77 —y oo ; 


where P is the classical Spalding heat transfer number If27l . This transfer number 
for slow vaporization belongs to the interval [0,0.1] and varies from P = (9(1) to 
P rs 20 for strong burning. This problem has been studied recently by Lu and Law 
[fl8l . These authors define an iterative method that has been shown to produce 
more accurate numerical results than the classical approximate solutions. 
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3.1 The non-ITM 


In the present case we consider the extended scaling group 

/* = A/, n* = \- 1 ri, P* = X- 2 P. (9) 


Let us notice that the governing differential equation and the two boundary condi¬ 
tions at r) = 0 in ([8]) are left invariant under ([9]) and, on the contrary, the asymptotic 
boundary condition is not invariant. By setting a value of P*, we can integrate the 
Blasius governing differential equation in ([8]) in the star variables on [0, T] r *j with 
initial conditions 


/*«>) = -/>*, j^(°) = 0 . |£( 0 ) = i, < 10 > 

in order to compute (r]^) & ^*(°°). Here J]£ is a suitable truncated boundary. 
The value of A can be found by 


A = 


dP 

dp 


(V*) 


1 1/2 


(ID 


After using (fill) to get the value of A, we can apply the scaling invariance to obtain 
the missing initial conditions 


/(0) = A-V, 0(0) -A’ 3 . (12) 

For the reader convenience, in table [2] we list sample numerical results. The 
case P* — P = 0 is, again, the Blasius problem d8]). In this case our non-ITM 
reduces to the original method defined by Topfer [j5J. The obtained skin friction 
coefficient is in good agreement with the values available in literature, see for 
instance the value 0.46599988361 computed by Fazio Il22ll . On the other hand, 
our value is different from the value 0.490 obtained by a 3 — 2 iteration solution 
of Lu and Law tfT8l . For the numerical results reported here, depending on the 
behaviour of the numerical solution, we have used 77 ^ = 5 or rj^ — 10. 

Figure [3] shows a sample numerical integration for P* — 1 that is transformed 
under © to P « 3.726397. We notice that the solution of the Blasius problem 
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p * 

dr H 

dr]* K ’ 

-/(0) 


P 

0 . 

1.655301 

0 . 

0.469553 

0 . 

0.1 

1.793644 

0.074668 

0.416289 

0.179364 

0.25 

2.025902 

0.175643 

0.346795 

0.506476 

0.5 

2.485809 

0.317129 

0.255152 

1.242904 

0.75 

3.048481 

0.429556 

0.187877 

2.286361 

1 . 

3.726397 

0.518031 

0.139016 

3.726397 

1.25 

4.528469 

0.587401 

0.103770 

5.660586 

1.5 

5.469166 

0.641403 

0.078184 

8.203749 

1.75 

6.548781 

0.683845 

0.059670 

11.460366 

2 . 

7.779561 

0.717055 

0.046086 

15.559122 

2.2 

8.863956 

0.738939 

0.037893 

19.500704 


Table 2: Surface gasification boundary condition: non-ITM results. 

with surface gasification boundary condition is computed by rescaling. Moreover, 
by rescaling we get r]^ < 7]^. 

In figure |4] we plot /(0) and j4(0) versus P. We know, from the literature, 
that as P goes to infinity then /(0) goes to —0.876. Moreover, as it is easily seen, 

j 2 r 

as P goes to infinity then j4(0) goes to zero. 

4 Concluding remarks 

The main contribution of this paper is the extension of the non-ITM proposed 
by Topfer |Q for the numerical solution of the celebrated Blasius problem 0 to 
classes of problems depending on a parameter. By requiring the invariance of the 
involved parameter we are able to solve the given BVP non-iteratively but for a 
different value of the parameter. This kind of extension was considered first by 
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Figure 3: Surface gasification boundary conditions with P* — 1. Numerical solu¬ 
tion by the non-ITM. 


Na IfTTfl . see also NA |fT2l Chapters 8-9]. Here we defined a non-ITM for Blasius 
equation with moving wall or surface gasification. 

Finally, the reader should be advised that non-ITM cannot be applied to all 
problems of boundary layer theory. In fact, sometimes we have to face problems 
that are not invariant with respect to all scaling groups. As an example we can 
consider the Falkner-Skan model lf28l 
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where / and 77 are similarity variables and P is a parameter related to the func- 
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Figure 4: Blasius problem with surface gasification boundary conditions. Non- 
ITM: behaviour of /(0) and j4(0) with respect to P. 

tional form od the fluid mainstream velocity. If we test the invariance of the gov¬ 
erning differential equation in (fT3l) under the extended scaling group 

T 7 * = A a, ?7 , /* = , P* = X a3 P, (14) 

where A is, again, the group parameter and 05/, for j — 1,2,3, are constant to be 
determined, then we get three invariant conditions 

a 2 -3ai =2(d5 2 -Q5i) = 05 3 = oj 3 + 2(a 2 - oti). (15) 

Now, it is a simple matter to show that the linear system defined by (fl5l) has 
the unique solution cq = a 2 = 05 3 = 0. However, an iterative extension of our 
transformation method has been developed in [f29l 30] and successfully applied to 
the Falkner-Skan model I29.I3H. 
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